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摘 要 : 基于 卫星 遥感 数据 ,提取 近 60 a 伊 塞 克 湖 面积 .水 位 变化 信息 , 反 演 伊 塞 克 湖 水 量变 化 时 间 序 列 ,结合 
1960 一 2020 年 CRU 和 气象 数据 .1960 一 2000 年 乔 尔 莲 阿 塔 气象 站 气温 降水 观测 数据 和 入 湖水 量 观测 数据 ,建立 湖泊 


水 量 平衡 模型 ,分析 水 量 平 衡 各 分 量 的 变化 特征 ,并 探讨 其 影响 因素 。 结 果 表 明 :(1) 1960 年 以 来 伊 塞 克 湖 水 量变 


化 经 历 了 持续 减少 -波动 增加 的 过 程 ,1998 年 为 变化 的 时 间 拐 点 ;20 世 纪 60 一 80 年 代 中 期 ,入 湖水 量 主要 受 灌溉 引 


水 影响 持续 减少 ,1986 年 后 随 灌溉 水 量 减少 .降水 和 冰 


中 融 水 的 增加 而 转 为 上 升 趋势 ;湖区 降水 以 9.1 mm- (10a) ”的 


速率 增加 ,蒸发 量 随 湖 区 升温 和 湖 体面 积 增加 总 体 呈 显著 增加 趋势 .2) 20 世 纪 80 年 代 中 期 以 前 伊 塞 友 湖 大 部 分 


年 份 湖泊 水 量 呈 负 平 衡 , 地 下 水 持续 补给 湖泊 ,1986 年 起 湖泊 的 水 量 收 支 亏损 逐渐 减 小 ,1998 年 以 来 以 正平 衡 为 
主 -(3) 入 湖 径 流 、 降 水 .蒸发 等 水 量 平 衡 分 量 的 互动 关系 决定 了 湖泊 水 量 的 变化 ,而 产 流 区 气候 变化 和 灌区 灌溉 引 
水 通过 改变 入 湖 径流 间接 驱动 湖泊 水 量 的 变化 ;1960 一 1986 年 ,以 灌溉 引水 为 主 的 人 类 活动 是 驱动 伊 塞 克 湖 水 量 


变化 的 主导 因素 ,贡献 率 达 71.6% ,1987 年 以 来 ,气候 变化 因子 对 湖泊 水 量变 化 的 累计 贡献 超过 80%。 
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亚洲 中 部 干旱 区 是 全 球 内 陆 湖泊 分 布 密集 的 
地 区 之 一 ,其 中 1 km 以 上 的 湖泊 有 3000 多 个 ,100 
km 以 上 的 大 型 湖泊 60 多 个 拉 。 在 全 球 气候 变 暖 的 
大 背景 下 ,中 亚 气 候 变化 及 其 水 文 响应 研究 得 到 了 
科学 界 的 广泛 关注 中 。 特 别 是 中 亚 干旱 区 的 内 陆 湖 
泊 ,不仅 是 区 域 气候 变化 的 指示 器 ,也 是 水 循环 的 
重要 组 成 部 分 ,探讨 内 陆 湖泊 水 量 收 支 变化 及 其 对 
变化 环境 的 响应 ,对 认识 干旱 区 气候 水 文 变化 规律 
有 重要 意义 。 

近 几 十 年 来 ,中 亚 和 干旱 区 内 陆 湖泊 发 生 了 剧烈 
的 变化 ,并 存在 明显 差异 2 ,而 影响 内 陆 湖泊 变化 的 
原因 较为 复杂 , 既 有 气候 波动 的 因素 ,也 有 人 类 活 
动 的 因素 ,同时 还 受 湖泊 大 小 ,形态 ,海拔 高 度 、 下 
垫 面 特征 和 水 文 循 环 等 影响 。 总 的 来 说 ,高 山 封 闭 
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湖 位 于 高 山 的 低洼 盆地 之 中 ,部 分 湖泊 具有 稳定 的 
冰川 融 水 补给 , 受 人 类 活动 干扰 较 小 ,对 区 域 气候 
变化 响应 敏感 ,如 赛 里 木 湖 . 阿 牙 克 库 木 湖 和 阿 其 
克 库 勒 湖 等 ;而 在 平原 区 ,内 陆 河 流域 农业 灌溉 引 
水 和 蓄 水 改变 了 水 资源 的 时 空 分 布 , 使 得 原来 可 以 
流入 下 游 尾 间 湖泊 的 部 分 或 全 部 径流 消耗 在 了 平 
原 绿洲 区 ,导致 旱 区 尾 闻 湖泊 水 量 平衡 发 生 剧 烈 变 
化 , 当 人 类 活动 的 干扰 大 于 气候 因素 的 影响 时 , 则 
表现 为 面积 萎缩 水 位 降低 和 水 量 减 小 ,如 威海 、 博 
斯 腾 湖 和 阿拉 湖 等 。 伊 塞 殉 湖 是 中 亚 最 大 的 高 山 
封闭 湖泊 ,发 源 于 昆 格 阿拉 套 山 和 杰 斯 克 耶 山 的 
100 多 条 河流 流向 伊 塞 克 湖 ;19 世 纪 后 期 以 来 , 随 着 
引水 垦荒 ,许多 河水 在 未 流入 湖 之 前 被 用 来 灌溉 田 
地 。 伊 塞 克 湖 的 变化 体现 了 气候 变化 因素 和 人 类 
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活动 干扰 的 复合 影响 , 且 具 有 明显 的 阶段 性 特征 。 
已 有 研究 表明 ,近代 以 来 伊 塞 克 湖 水 位 总 体 呈 下 降 
趋势 一 ,变化 过 程 可 以 分 为 波动 下 降 - 快 速 下 降 - 组 
慢 波动 上 升 3 个 阶段 ;湖泊 面积 变化 与 水 位 变化 
有 较 强 的 同步 性 , 呈 先 减 小 后 增加 的 趋势 "2 20 
世纪 80 年 代 以 前 引水 灌溉 导致 了 伊 塞 克 湖 水 位 持 
续 下 降 “ ”1987 年 以 来 ,中 亚 地 区 气候 发 生变 
化 ,主要 表现 为 降水 量 .冰川 融 水 量 和 河川 径流 量 显 
著 增 加 ,间接 促进 了 湖泊 水 位 上 升 汪 。 然 而 ,以 往 研 
究 大 多 分 析 面 积 或 水 位 单一 水 文 要 素 变 化 及 其 驱 
动因 素 %", 且 数据 主要 集中 于 2010 年 以 前 ,缺少 
对 伊 塞 克 湖 近期 变化 的 研究 ,特别 是 对 湖泊 水 量 平 
衡 分 量 和 收 支 状态 的 变化 ,以 及 水 量 平衡 各 分 量 对 
湖泊 水 量 的 定量 贡献 等 尚 不 明确 。 

因此 ,本 文选 择 伊 塞 克 湖 作 为 研究 区 ,基于 
1990—2020 年 Landsat TM/ETM +/OLI , Jason1/2/3 等 
遥感 影像 数据 ,提取 湖泊 面积 水 位 信息 ,结合 历史 
观测 数据 ,重建 长 时 序 面积 和 水 位 变化 序列 , 反 演 
伊 塞 克 湖 水 量 ; 选 取 伊 塞 克 湖 流域 乔 尔 途 阿 塔 气象 
站 1960 一 2000 年 气象 数据 对 1960 一 2020 年 CRU 再 
分 析 气 象 数 据 进行 修正 ,建立 湖泊 水 量 平衡 模型 ， 
定量 分 析 降 水 .蒸发 、 人 湖 径 流 与 湖泊 水 量 等 多 要 
素 变 化 及 其 互动 关系 ,分 析 气 候 变化 和 人 类 活动 关 
键 因子 对 伊 塞 克 湖 水 量变 化 的 贡献 。 研 究 结果 可 
为 伊 塞 克 湖 流域 水 资源 调控 和 生态 保护 提供 数据 
MTF, 


1 研究 区 概况 


伊 塞 克 湖 位 于 吉尔 癌 斯 斯 坦 天 山北 坡 的 封闭 性 
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山 间 伟 地 中 (图 1), 属 内 陆 咸水湖 。 湖 心 位 于 77.33°E， 
42.42°N 附近 ,湖面 平均 海拔 1607.03 m ,湖岸 线 周 长 
约 669 km ,平均 深度 278.4 了 ,最 大 深度 约 702 m, WH 
水 体积 17350x10 mm ,总 面积 约 6236 km 。 湖 区 
位 于 大 陆 性 气候 带 中 部 ,气候 温和 干燥 ,1 月 平均 气 
温 -6% ,7 月 平均 气温 15~25 % ,湖面 年 均 蒸 发 量 约 
820 mm, 年 降水 200~300 mm。 湖 周 山谷 环抱 , 北 侧 
为 Kungey Alataw 山脉 ,海拔 4700 m; 南 侧 为 Terskey 
Alataw 山脉 ,最 高 海拔 5200 m。 这 些 山脉 围 成 了 约 
22080 km 的 集 水 流域 "。 流 域内 冰川 面积 约 650.4 
km ,有 102 条 河流 汇 人 该 湖 。 伊 塞 克 湖 为 高 山 封 闭 
内 陆 湖 ,许多 河水 在 未 汇 人 湖 之 前 就 被 用 于 灌溉 等 
生产 活动 ,但 河流 不 是 唯一 的 补给 来 源 ,每 年 湖泊 
平均 地 表 水 补给 量 约 为 13 x 10° ms ,平均 地 下 水 补 
给 量 约 为 20x10 m, 


2 材料 与 方法 


2.1 数据 收集 

基于 30m 分 辩 率 的 Landsat TM/ETM +/OLI 2% 
影像 提取 了 1990 一 2020 年 伊 塞 克 湖面 积 变 化 时 间 
序列 ,影像 时 相 集 中 在 5 一 8 月。 湖泊 水 位 数据 的 获 
取 分 为 2 个 阶段 :1960 一 1991 年 和 1992 一 2020 年 。 
前 一 阶段 数据 来 自 Romanovsky 的 研究 ,后 一 阶段 
水 位 信息 从 TOPEX/Poseidon(T/P) 卫 星 以 及 Jason1/2/ 
3 测 高 卫星 数据 产品 中 提取 *。 

选取 距离 伊 塞 克 湖 最 近 的 气象 站 乔 尔 鞍 阿 塔 
(海拔 1645 mm) 为 代表 站 ,气象 站 点 数据 来 自 国家 
冰川 冻 土 沙漠 科学 数据 中 心 (http:/www.ncdc.ac. 
cn) Y ,数据 参量 为 逐 月 平均 气温 和 逐 月 降雨 量 ， 
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图 1 研究 区 概况 


Fig. 1 Overview of the study area 
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时 段 为 1960 一 2000 年 。2001 一 2012 年 月 平均 气温 
和 月 降水 数据 由 伊 丽 努 尔 阿 力 甫 江 等 "的 研究 提 
供 。 因 地 理 条 件 限制 及 政策 原因 ,相关 数据 缺失 严 
重 ,需要 利用 气象 再 分 析 数 据 进 行 补 全 。 采 用 英 
East Anglia 大 学 Climatic Research Unit(CRU) 发 布 的 
CRU TS 4.05 数据 (http://data.ceda.ac.uk/badc/cru/da- 
ta/cru_ts/) , 23 BIIEN 0.5°x0.5° ,时 段 为 1960 一 
2020 4 , ,所 用 数据 参量 包括 逐 月 平均 气温 、 逐 月 降 
水 量 。 

土地 覆盖 数据 来 自 欧洲 太空 局 Climate Change 
Initiative (CCI) (https://www.esa-landcover-cci.org) , 空 
间 分 辩 率 为 300 m x 300 m, 用 于 提取 1990—2015 年 
伊 塞 克 湖 流域 耕地 面积 变化 时 序 ;1960 一 2000 年 人 
湖 流量 观测 数据 由 王国 亚 等 的 研究 提供 ;2001 一 
2017 年 径流 数据 来 自 Uwamahoro Se XY} FASE Fe WH 
流域 水 文 模拟 的 结果 ;湖面 蒸发 观测 数据 来 自 
CAwater (http://www.cawater-info.net/index_e.htm ) ; YK 
川 面积 数据 来 自 中 国 科 学 院 “ 地 球 大 数据 科学 工程 
数据 共享 服务 系统 ”提供 的 1990 年 .2000 年 .2010 
年 .2015 年 4 期 中 亚 土地 利用 与 覆盖 数据 (http:/da- 
ta.casearth.cn/) 2 ;研究 区 水 系数 据 来 自 国家 冰川 冻 
土 沙 漠 科 学 数据 中 心 ”; 伊 塞 殉 湖 流域 年 平均 灌溉 
耗 水 量 来 自 王国 亚 等 ”的 研究 ,数据 时 段 为 1935 一 
2000 年 。 
2.2 湖泊 水 量 计算 

基于 1990 一 2020 年 Landsat TM/ETM+/OLI 1/8 
影像 ,采用 “全 域 - 局 部 ” 自 适应 迭代 的 国 值 分 割 方 
法 提取 湖泊 面积 信息 ” ,并 加 以 人 工 修订 ,得 到 
1990 一 2020 年 多 期 湖泊 面积 变化 序列 。 由 于 伊 
塞 克 湖面 积 与 水 位 变化 的 同步 性 较 强 ,可 根据 
Romanovsky'“ 的 研究 提供 的 历史 水 位 数据 , 推 求 
出 1960 一 1989 年 间 的 面积 变化 时 间 序 列 。 根 据 
遥感 提取 的 水 位 .面积 时 间 序 列 ,估算 伊 塞 克 湖 
REKE 。 
2.3 湖面 蒸发 计算 

采用 位 于 伊 塞 克 湖 北 侧 的 乔 尔 莲 阿 塔 站 实测 
数据 (气温 .日照 时 数 等 ) 计 算 湖面 蒸发 。Penman- 
Monteith 公式 无 论 在 湿润 区 还 是 在 干旱 区 表现 较 
好 ” ,但 所 需 参 数 包 括 最 高 气温 .最低 气温 、 风 速 、 
日 照 时 数 .相对 湿度 等 ,对 气象 资料 的 完整 性 有 很 
高 的 要 求 , 一 定 程度 上 限制 了 该 公式 的 使 用 ;Harg- 
reaves 方法 仅 基于 温度 和 地 理 位 置 即 可 计算 蒸 散 


发 , 且 经 过 校正 后 结果 较为 合理 ”。 本 研究 使 用 
Hargreaves 公 式 ( 以 下 简称 也 氏 公 式 ) 计 算 湖 面 蒸 发 
量 ,H 氏 公式 所 需 参 数 有 日 平均 最 高 气温 日 平均 最 
低 气 温 以 及 外 空 辐射 ”, 其 中 日 平均 最 高 气温 和 日 
平均 最 低 气 温 数 据 由 气象 站 提供 ,外 空 辐射 资料 可 
以 根据 纬度 查阅 FAO 相关 文献 得 到 。 
W=A x ET 


(1) 
ET = 0.0023R (T+ 17.8)7 Tya, (2) 
式 中 :下 为 湖泊 蔡 发 总 量 (m ); 4 为 湖泊 的 面积 
(m); ET 为 应 用 HH 氏 公式 计算 的 参考 蒸 散 量 
(mm-d'); R, 为 与 外 空 辐射 等 价 的 水 量 (mm* d); T 
为 日 平均 气温 (% ) ; 7 为 日 平均 最 高 气温 (%C); 
7T ， 为 日 平均 最 低 气 温 (% ) ;0.0023 为 初始 公式 的 经 
验 值 。 
2.4 降水 和 气温 时 间 序 列 重 建 
站 点 实测 数据 精度 高 ,但 在 反映 气温 、 降 水 时 
空 特征 方面 有 局 限 性 ;CRU 格 网 数据 能 够 反映 降水 
和 气温 的 时 空 分 布 特征 ,但 是 有 一 定 的 系统 误差 。 
此 有 必要 对 格 点 气温 和 降水 数据 进行 修正 ,弥补 
FE”, 
本 文 将 1960 一 2000 年 作为 拟 合 期 ,2001 一 2012 
年 作为 验证 期 ,采用 逐 月 比例 缩放 法 "来 修正 
2013 一 2020 年 时 段 内 CRU 降 水量 序列 ,公式 如 下 : 
a,P,=0, (3) 
Poor = a(p;) (4) 
式 中 世 为 月 份 , 取 1,2，……,12; a, 、 5 分 别 为 实测 与 
模拟 的 月 降水 比值 和 平均 比值 ; P, 为 拟 合 期 CRU 
月 平均 降水 量 , mm ; O, 为 拟 合 期 气象 站 点 月 平均 
降水 量 , mm ; 疡 为 未 修正 时 段 CRU 月 降水 量 ， 
mm ; 五, 为 修正 后 的 CRU 月 降水 量 , mm o 
为 了 评价 误差 校正 后 的 再 分 析 数 据 在 研究 区 
伊 塞 克 湖 流域 的 适用 性 ,采用 相对 误差 指标 和 Pear- 
son 相关 系数 ( 尼 ) 作 为 评判 标准 。 相 对 误差 (ARE ) 
计算 公式 为 : 
Non -N 
ARE - w 


obs 


x 100% (5) 


式 中 : NV,,,, 为 模拟 值 ; N。 为 观测 值 。 

温度 序列 重建 采用 线性 比例 缩放 法 ,将 
1960 一 2000 年 时 段 内 站 点 月 平均 气温 作为 实测 值 ， 
将 CRU 月 平均 气温 作为 模拟 值 ,公式 如 下 : 


= T+ (Mins fiz.) (6) 
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式 中 : 7, 为 校正 后 气温 , C; 7 为 校正 前 模拟 气 
温 , Cos jw 为 站 点 实测 气温 , C ; ,为 模拟 气温 
值 , Co 

将 2001 一 2012 年 时 段 内 站 点 数据 与 校正 后 的 
CRU 数据 对 比分 析 ,检验 校正 后 的 CRU 数据 的 有 
效 性 。 

从 表 1 可 以 看 出 ,经 过 校正 的 CRU 降水 数据 与 
站 点 数据 的 相关 性 有 所 提高 , 记 从 0.63 提高 到 0.74， 
相对 误差 从 30.91% 下 降 到 12.51% ,经 过 校正 后 的 
CRU 温度 数据 与 实测 年 均 气 温 误差 在 5% 以 内 ,说 
明 经 过 校正 后 的 CRU 降水 和 气温 序列 能 达到 精度 
HORO ,校正 结果 可 信和 度 高 ,该 重建 时 间 序 列 方法 
在 研究 区 适用 。 


R1 温度 和 降水 时 序 重 建 结果 精度 评价 
Tab.1 Accuracy evaluation of temperature and 


precipitation time series reconstruction results 


校正 前 校正 后 
R ARE/% R ARE/% 
气温 0.78 104.86 0.78 3.71 
降水 0.63 30.91 0.74 12.51 


25 水 量 平衡 模型 
某 时 段 湖 泊 库 容 差 与 增 ( 减 ) 水 的 关系 ,可 用 水 
量 平衡 方程 式 *” 表示: 


AV _ 7 _ 
AV =A(h\(P-E)+04- Om D 


式 中 : A 为 计算 时 段 (a); AV 为 湖泊 库容 变化 量 
(m°); 4 为 湖泊 水 面 面 积 (mw), A 是 水 位 h AY PRR 
P 为 研究 区 内 降雨 量 ( mm ) ;为 研究 区 内 总 共 发 
量 ( mm ); Q,, HAWKE Mm); Qua 为 出 湖水 量 


(mê) 


1609.0 F Ca) 一 水 位 — 面积 6230 
1608.5 6220 


1608.0 


1607.5 


水 位 mm 


1607.0 


1606.5 


À 
1606.0 4 ! L 1 L 1 1! 6180 


1960 1970 1980 1990 2000 2010 2020 
年 份 
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由 于 伊 塞 克 湖 属于 封闭 湖泊 ,没有 出 水 河流 ， 
根据 伊 塞 克 湖 流域 实际 情况 ,将 水 量 平衡 方程 表 
INN: 

WV = A(h)(P-£)+Q 

=A(h)(P-E) + Quia + Q sess 
式 中 : Qaira 为 河流 入 湖水 量 (m); Qara + 二 下 出 
为 地 下 水 入 .出 湖水 量 (m )。 

余 项 表示 为 : 

Qan = Q wra T Qiu ri (9) 


3 结果 与 分 析 


3.1 伊 塞 克 湖面 积 、 水 位 与 水 量变 化 

伊 塞 克 湖 面积 与 水 位 变化 的 同步 性 较 强 , 呈 先 
下 降 后 上 升 的 总 趋势 (图 2)。 如 图 2a 所 示 ,1960 一 
1986 年 , 伊 塞 克 湖面 积 快速 减 小 ,1987 年 起 湖泊 面 
积 减 小 的 趋势 减缓 ,但 由 于 1997 年 流域 降水 量 为 近 
60 a 最 低 , 仅 为 148.1 mm( 图 3) ,导致 1998 年 湖泊 面 
积 较 1960 年 减少 34.91 km ,水 位 降 至 1606.17 m; 
1998 年 后 随 着 降水 量 增 加 湖泊 水 位 和 面积 变化 转 
为 上 升 趋势 ,至 2011 年 面积 持续 增加 至 6213.72 
km2? ,水 位 增加 到 1607.84m ,此 时 湖泊 水 位 基本 恢复 
到 1960 年 前 后 水 平 ;2011 一 2020 年 湖泊 面积 水 位 
小 幅 波动 下 降 , 维 持 在 6210.09 km? 和 1607.63 m 
左右 。 

伊 塞 克 湖 属 于 高 山内 陆 湖 泊 , 湖 盆 由 地 壳 陷 落 
形成 , 横 形 湖 倪 决定 了 湖泊 水 量 的 变化 不 会 太 剧 
烈 。 如 图 2b 所 示 ,1960 年 以 来 ,湖泊 水 量 最 大 变化 
不 超过 86x10* m, 4HE F 16000 x10 wm 以 上 的 总 水 
量 来 说 ,变化 很 小 。 水 量 的 变化 趋势 与 水 位 .面积 


河流 入 十 Owen T Quen (8) 


16180 r (b) 
16160 | a 
16140 + w Pi 


16120 上 Me 


面积 /km 
水 量 /10m3 


a Py 

16100 上 Met 
ms = 

16080 上 a 


16060 1 1 1 1 1 1 jj 
1960 1970 1980 1990 2000 2010 2020 
年 份 


图 2 1960 一 2020 年 伊 塞 克 湖 水 位 .面积 和 水 量变 化 


Fig. 2 Changes in water level, area and volume of Issyk-Kul Lake from 1960 to 2020 
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的 变化 一 致 ,1960 一 1998 年 期 间 持 续 减少 ,减少 速 
率 约 2.1x10 ma!。1998 一 2006 年 间 入 湖 流量 每 
年 约 增加 0.084x10* m’, 湖 泊 水 量 处 于 缓慢 上 升 阶 
Be, 至 2006 年 增加 到 16105.7x10* m°, 2006—2008 
年 水 量 小 幅 减 少 9.4x10’ m,2009 一 2011 年 年 均 降 
水 量 334.26 mm, 高 于 多 年 平均 降水 量 约 54 mm, 7K 
量 迅 速 增加 ,2011 年 达到 16143.6x10° m°, 与 1966 
年 水 量 持 平 。 此 后 水 量 呈 小 幅 波 动 变化 ,变化 率 
为 -3.8x108 m°- (10a)'。 

3.2 伊 塞 克 湖 水 量 平衡 分 量变 化 

3.2.1 降水 变化 ， 伊 塞 克 湖 湖区 多 年 平均 降水 量 为 
279.8 mm; 1960 年 至 今 降 水 量 呈 增加 趋势 ,增幅 率 
9.1 mm-(10a)'(13). FEP 1993—1997 年 降水 量 显 
著 下 降 ,由 1993 年 的 420.2 mm 降低 到 1997 年 的 
148.1 mm ,也 是 近 60 a 降水 量 最 低 值 。1998 一 2005 
年 降水 量 以 3.7 mma 1! 的 幅度 持续 增加 。 此 后 呈现 
波动 增加 趋势 , 至 2020 年 降水 量 达 到 了 321.8 mm, 
是 多 年 平均 降水 量 的 1.2 倍 左右 。 


450 一 降水 量 — 线性 
y=0.9058x-1522.6 = 


降水 量 /mm 


1960 1970 1980 1990 


年 份 


2000 2010 2020 


图 3 1960 一 2020 年 伊 塞 克 湖 降水 量变 化 
Fig. 3 Precipitation change in Issyk-Kul Lake from 
1960 to 2020 


3.2.2 湖面 蒸发 变化 HAA WAR W 
区 气温 变化 有 较 强 的 同步 性 。1960 一 2017 年 间 ,多 
年 平均 蒸发 量 约 927.6 mm, 与 CAwater 发 布 的 观测 
蒸发 量 相近 ;湖面 蒸发 量 总 体 呈 显著 增加 趋势 ( 显 
著 性 水 平 a=0.01) ,增幅 为 15.7 mm-(10ay'. 1960— 
1985 年 ,湖泊 水 面 蒸 发 维持 在 较为 稳定 的 状态 , 平 
均 约 911.7 mm; 1986 一 1993 年 间 蒸 发 量 波动 较 大 ; 
1993 年 后 随 着 湖区 气温 升 高 ,湖面 蒸发 波动 增加 ; 
1994—2002 4F4F IY AR IA 972.5 mm;2003 一 2017 


a 39 卷 
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À - 1200 
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图 4 1960 一 2017 年 伊 塞 克 湖 蒸发 量变 化 
Fig. 4 Changes in evaporation from Issyk-Kul Lake, 
1960-2017 


AE AR A AV 5.8 mma” 的 速率 显著 增加 (显著 性 水 
平 a=0.01 )。 

3.2.3 入 湖 径流 如 网 5 所 示 , 近 60a 伊 塞 克 湖人 湖 
流量 总 体 呈 增加 趋势 ,年 际 变化 范围 在 30.23x108 
m°~49.62x 10° m? Z [i] ,多 年 平均 入 湖 径流 量 为 
39.41x10'"m。 由 M-K 检 验 结果 可 知 ,1986 年 人 湖 径 
流量 发 生 了 由 减少 到 增加 的 突变 ,这 一 时 间 拐 点 与 
中 亚 气 候 由 “ 暧 干 ” 向 “ 暧 湿 ” 转 变 的 突变 时 间 基 本 
一 致 和。1986 年 流域 降水 量 迅速 增加 , 相 较 1985 年 
增加 了 6.27x10: ms; 加 上 升温 导致 冰雪 融 水 对 径流 
的 补给 增加 ,是 入 湖 流量 由 持续 减少 转变 为 增加 趋 
势 的 主要 原因 。 从 2 个 阶段 来 看 ,1960 一 1986 年 伊 
塞 克 湖 流域 降水 量 总 体 处 于 较 低 水 平 , 和 湖 流量 主 
要 受 农 业 灌溉 引水 影响 持续 减少 。 王 亚 国 等 “” 指 
出 ,1950 一 1980 年 代 中 期 灌溉 用 水 量 在 伊 塞 克 湖 盆 
地 耗 水 量 中 占 比 最 大 ,扰乱 了 地 下 水 和 地 表 水 的 自 
然 交 换 。1987 年 起 ,流域 降水 在 大 部 分 年 份 维持 在 
较 高 水 平 ;1989 一 1992 年 农业 .工业 ,城市 用 水 量 分 
别 下 降 了 66% . 36% .64%55 ,流域 耗 水 总 量 与 
1977—1989 年 相 比 减少 了 53.49% ; 直至 2015 年 研 
究 区 内 灌溉 耗 水 量 处 于 稳定 状态 。 径 流 补给 量 激 
增 与 灌溉 引水 减少 的 双重 作用 使 得 入 湖 流量 维持 
了 逐步 回升 的 趋势 ,至 2008 年 入 湖 流量 达到 了 
49.07x10* m°; 2010—2015 年 降水 量 减少 ,入 湖 流量 
也 随 之 下 降 了 12.8x10* ms ,此 后 受降 水 回升 和 冰川 
消融 的 影响 ,入 湖 径流 量 迅速 增加 ,到 2017 年 人 湖 
径流 量 达 到 44.61x10: m ,高 于 多 年 平均 入 湖 径 流量 
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5 1960 一 2017 年 伊 塞 克 湖 入 湖 流量 变化 
Fig. 5 Change of the inflow to Issyk-Kul Lake from 1960 to 2017 


443 5.2x10° m*o 

3.2.4 余 项 依据 水 量 平衡 公式 , 余 项 可 近似 为 地 
下 水 的 变化 量 ,反映 伊 塞 克 湖 与 地 下 水 的 互动 关 
系 。 如 图 6a 所 示 ,1960 一 1986 年 间 大 部 分 年 份 伊 塞 
克 湖 水 量 平衡 余 项 为 负 值 ,湖泊 显示 了 负 平 衡 , 表 
明 地 下 水 持续 补给 湖泊 , 除 1964 年 、1969 年、1982 
年 外 ,地 下 水 对 湖泊 的 平均 补给 量 约 4.5x10: m-a, 
其 中 1976 年 地 下 水 补给 最 多 ,补给 量 13.2x10’ m°; 
除 1990 年 外 ,1986 一 1994 年 间 水 量 平衡 余 项 均 为 正 
值 ,湖泊 水 量 收 支 为 正 ,表明 这 些 年 份 以 湖泊 补给 
地 下 水 为 主 。1995 一 2020 年 , 余 项 在 -12.7x10* mi~ 


水 量 平衡 余 项 /108ms 


13.2x10° m 之 间 波 动 变化 。 

伊 塞 克 湖水 量 平衡 余 项 的 年 代 变 化 如 图 6b 所 
示 ,20 志 纪 80 年 代 以 前 ,湖泊 水 量 基本 处 于 负 平 衡 
状态 ,地 下 水 持续 补给 湖泊 ,其 中 70 年 代 地 下 水 补 
给 较 多 ,平均 为 4.76x10 m a';20 世 纪 80 年 代 到 21 
世纪 10 年 代 湖 泊 水 量 处 于 三 余 的 状态 ,地 表 径 流 和 
降水 的 补给 能 够 维持 湖泊 的 水 量 平衡 。 

由 图 7 可 知 ,1980 年 以 前 降水 .入 湖 流量 .蒸发 
量 和 水 量 平衡 余 项 均 为 负 距 平 ,说 明 该 时 段 水 量 平 
衡 各 分 量 均 偏 低 ,湖泊 水 量 处 于 负 平 衡 状 态 ,湖泊 
水 位 持续 下 降 ;1980 一 1989 年 降水 量 为 正 距 平 ,入 
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图 6 1960 一 2017 年 伊 塞 克 湖水 量 平 衡 年 际 变化 (a) 和 年 代 变化 (b) 
Fig.6 Inter-annual (a) and chronological (b) change of water balance in Issyk-Kul Lake, 1960-2017 
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5ST 四 降水 SS 人 湖 流量 co Roki m 水 量 余 项 倾向 率 达 0.369 CC.(10a) ,变化 幅度 (-1.3~1.5 %C ) 远 
y “T a 大 于 全 球 气温 的 变化 幅度 (-0.6~0.55 C). YKJ 
è 3f 处 于 持续 的 负 物 质 平衡 而 不 断 退缩 ,亚洲 高 山区 负 
22 物质 平衡 (-150+110 kg.m*…a')""。 围 绕 伊 塞 克 湖 
=} 的 100 多 个 山区 子 流域 是 主要 产 流 区 ,流域 气候 变 
化 趋势 与 中 亚 地 区 的 整体 趋势 具有 较 强 的 一 致 
y 性 。1960 年 以 来 , 产 流 区 气温 呈 波 动 上 升 趋势 显 


图 7 1960 一 2017 年 伊 塞 克 湖水 量 平衡 各 分 量 年 代 距 平 
Fig.7 Chronological distance of each component of the 
water balance of Issyk-Kul from 1960-2017 


湖 流量 和 蔡 发 量 为 负 距 平 ,表明 期 间 降 水 增多 ,入 
湖 流量 减少 ,降水 量 的 增加 抵消 了 入 湖 流量 减少 的 
部 分 影响 ,加 上 湖面 莹 发 量 减 少 导 致 水 量 平衡 余 项 
为 正 距 平 ,期 间 湖泊 水 量 处 于 鱼 余 状态 ,对 应 1986 
年 后 伊 塞 克 湖 水 位 不 再 继续 下 降 , 而 是 在 波动 中 维 
持 稳 定 ( 图 2a);1990 一 2009 年 湖泊 水 量 平衡 各 分 量 
均 显 示 为 正 距 平 ,期 间 湖泊 水 位 回升 ;2010 一 2017 
年 降水 量 `. 和 湖水 量 和 蒸发 量 距 平 为 正 , 但 湖面 蒸 
发 的 增 量 超过 降水 量 与 人 湖 流量 的 增 量 之 和 ,导致 
这 一 阶段 湖泊 水 量 处 于 负 平 衡 状态 ,地 下 水 少量 补 
给 湖泊 。 


4 讨论 


伊 塞 克 湖 是 地 处 高 山 盆地 的 封闭 内 陆 湖 ,湖水 
来 源 包括 大 小 河流 的 汇 人 、 降 水 以 及 地 下 水 的 补给 
(水 量 收入 ) , 湖 体 水 面 蔡 发 是 水 量 耗 散 的 主要 途径 
(水 量 支出 ); 而 从 流域 尺度 来 看 , 产 流 区 的 气候 变 
化 和 湖 周 绿洲 区 以 灌溉 引水 为 主 的 人 类 活动 ,通过 
影响 入 湖 流量 间接 驱动 湖泊 水 量 的 变化 。 

沈 永 平等 的 研究 指出 ,北大 西洋 涛 动 与 欧 亚 
大 陆 北 部 上 空 的 西风 环流 正 相 关 , 北 大 西洋 涛 动 指 
数 1956 一 1971 年 是 持续 负 位 相 期 ,1972 一 1986 年 处 
于 正 负 波动 期 ,而 1987 一 2000 年 基本 是 持续 正 位 相 
期 ,水 汽 增加 可 能 是 中 亚 地 区 降水 增加 的 主要 原 
因 。 与 此 同时 ,20 世 纪 80 年 代 初 期 ,中 亚 气温 变化 


著 , 量 1986 年 后 升温 速率 明显 高 于 前 期 ;1960 一 
1985 年 间 降 水 量 呈 波动 变化 ,1986 年 起 迅速 增加 且 
波动 剧烈 ,1986 一 2020 年 降水 量 较 前 期 平均 增加 
29.34 mm -a` (图 8a)。 随 着 高 山区 升温 加 剧 ,冰雪 
融 水 对 径流 的 补给 量 增加 ,气温 对 湖泊 变化 的 正 效 
应 凸显 下 。 一 些 学 者 基于 典型 冰川 的 观测 数据 , 揭 
示 了 伊 塞 克 湖 流域 冰 川 物质 平衡 减少 的 事实 ”1。 
1958 一 2016 年 伊 塞 克 湖 北部 的 土 尤 克 苏 冰川 表面 平 
均 下 降 23.2+2.2 m, 消 失 了 约 67.7+6.7x105 m“, 
Romanovsky“ 指 出 1973 年 开始 ,位 于 伊 塞 克 湖 南 侧 
的 帖 尔 斯 克 伊 山 (Terskei Ala-Too) 上 的 卡拉 :巴特 卡 
克 (Kara-Bat-kak ) 冰 川 物质 平衡 为 负 , 这 期 间 所 有 基 
于 雪 、 冰 川 补给 的 河流 都 观测 到 了 冰川 径流 的 增 
加 , 约 4 km 冰川 体积 的 融化 增加 了 入 湖 径流 ,这 相 
当 于 流域 内 8.3% 的 冰川 。 本 文 提取 1990 年 .2000 
年 .2010 年 和 2015 年 伊 塞 克 湖 流域 冰川 面积 (图 
8b) ,结果 表明 1990—2015 年 流域 冰川 面积 持续 减 
少 , 由 1990 年 的 691 km 减少 至 2015 年 的 647 km’, 
减 小 速率 约 1.26 km…a'。 王 国 亚 等 中 的 研究 表明 ， 
1935—2000 年 伊 塞 克 湖 人 湖 径 流量 呈 增 加 趋势 , 变 
化 速率 为 9.2 mm: (10a)' ; 1975 年 以 后 增加 趋势 更 为 
明显 ,为 24.4 mm(10a)'。 这 与 1980 年 代 以 来 ,湖泊 
水 位 波动 回升 的 趋 热 相 一 致 ,说 明 产 流 区 降水 和 冰 
雪 融 水 增加 导致 入 湖 径 流量 增加 ,驱动 了 伊 塞 克 湖 
水 量 由 负 平 衡 到 正平 衡 的 转变 。 

以 灌溉 引水 为 主 的 人 类 活动 也 是 影响 伊 塞 克 
湖 入 湖水 量 的 重要 因素 。 伊 塞 克 湖 流域 土地 利用 
开发 可 追溯 至 19 世纪 后 期 至 20 世 纪 中 期 , 伊 塞 克 
湖 盆 地 引水 垦荒 2 。1930 年 流域 灌溉 面积 仅 为 
500 km2,20 世 纪 50 年 代 至 1980 年 代 中 期 灌溉 面积 
增 至 1540 km?) ; 这 与 1960 一 1986 年 期 间 和 人 湖 流量 
持续 减少 和 湖水 水 位 快速 下 降 相 一 致 。 施 雅 风 史 | 
指出 ,1935 一 1980 年 间 伊 塞 克 湖 地 表 和 地 下 径流 补 
给 明显 减少 ,灌溉 引水 量 .灌区 营 散 发 量 以 及 湖面 
蒸发 量 上 升 ,致使 伊 塞 克 湖 水 位 明显 下 降 。 邵 新 媛 ” 
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图 8 1960 一 2020 年 伊 塞 克 湖 流域 气温 、 降 水 量变 化 (a) 和 1990 一 2015 年 冰川 面积 变化 (b) 
Fig. 8 Changes in temperature and precipitation in the Issyk-Kul from 1960 to 2020 (a) and 


change in glacier area from 1990 to 2015 (b) 


的 分 析 结 果 表明 ,1927 一 1983 年 由 于 农田 用 水 .工业 
用 水 和 城市 生活 用 水 , 伊 塞 克 湖 人 湖水 量 每 年 减少 
约 6.0 10° mm, 导致 湖水 位 每 年 下 降 约 $ cm; 其 中 
1977—1989 年 间 灌 溉 耗 水 较 1935—1976 年 增加 
160%, 1989—1992 年 ,耕地 大 范围 扬 荒 ,农业 、 
工业 ,城市 用 水 量 分 别 下 降 了 66% 、36% .64% = 。 
本 文 基于 欧洲 太空 局 土地 利用 变化 数据 提取 
1992—2015 年 仇 塞 克 湖 流域 耕地 面积 变化 (图 9a)， 
结果 表明 :1992 年 以 来 仇 塞 克 湖 流域 耕地 面积 小 幅 
增加 了 184.56 km ,依据 伊 塞 克 湖 流域 净 灌 溉 定额 ”， 
结合 1990—2000 年 的 平均 耗 水 量 "" ,计算 得 到 
2001—2015 年 间 的 平均 耗 水 量 为 20.72x10* m° (图 
9b), 与 1990 一 2000 年 间 平 均 灌 溉 耗 水 量 基本 持平 。 

水 量 收入 与 支出 间 的 互动 关系 决定 了 伊 塞 克 
湖 的 历史 演变 进程 与 未 来 趋势 。 入 湖 径流 、 降 水 、 
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蒸 散 发 等 水 量 平 衡 分 量 的 互动 关系 决定 了 湖泊 水 
量 的 变化 ,而 产 流 区 气候 变化 和 灌区 灌溉 引水 通过 
改变 入 湖 径流 间接 驱动 湖泊 水 量 的 变化 。 因 此 , 选 
取 山 区 降水 和 气温 灌溉 引水 量 、 湖 区 降水 和 蔡 发 
为 主要 驱动 因子 ,通过 析 因 分 析 法 计算 各 因子 对 伊 
塞 克 湖水 量变 化 的 贡献 率 ( 表 2) ,结果 表明 :气候 变 
化 与 人 类 活动 对 伊 塞 克 湖 水 量变 化 的 影响 阶段 性 
差异 显著 ,1960 一 1986 年 以 灌溉 引水 为 主 的 人 类 活 
动 是 驱动 湖泊 水 量 急剧 减少 的 主导 因素 ,页 献 率 达 
71.6% ; 1986 年 以 来 ,气候 的 " 暖 湿 "化 变化 ,以 及 灌 
溉 引水 量 的 减少 ,促进 了 伊 塞 克 湖 水 位 持续 升 高 ， 
并 进入 稳定 恢复 的 状态 ,气候 变化 因子 的 累计 贡献 
达到 80% 以 上 。 值 得 注意 的 是 ,在 流域 持续 升温 ， 
且 湖 区 升温 幅度 较 大 的 情况 下 ,湖泊 面积 增加 必然 
导致 蒸发 量 不 断 增 大 ,湖区 莹 发 对 水 量变 化 的 影响 
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图 9 1992 一 2015 年 伊 塞 克 湖 流域 耕地 面积 变化 0) 和 1960 一 2015 年 伊 塞 克 湖 流域 灌溉 耗 水 量变 化 (b) 
Fig. 9 Changes in cropland area in Issyk-Kul basin from 1992-2015 (a) and changes in irrigation water consumption in Issyk-Kul 
basin from 1960-2015 (b) 
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表 2 各 因子 对 伊 塞 克 湖 水 量变 化 的 贡献 率 
Tab. 2 Contribution of each factor to the variation of 


water volume in Issyk-Kul Lake 


影响 因子 Z 

1960—1986 4E 1987—2017 4E 
山区 降水 15.17 22.11 
山区 气温 2.31 3.08 
灌溉 引水 71.64 13.85 
湖区 降水 10.34 15.72 
湖区 蒸发 0.55 45.25 


加 剧 。 这 种 情况 在 2010 年 后 表现 更 为 突出 ,蒸发 耗 
散 对 湖泊 水 量 的 钠 效 应 明显 超过 入 湖 径 流 和 降水 
的 正 效应 ,导致 了 近期 伊 塞 克 塞 克 湖面 积 ` 水 位 和 
水 量 的 小 幅 波 动 减 小 。 


5 结论 


本 文 以 伊 塞 克 湖 为 研究 对 象 ,基于 多 源 遥 感 数 
据 重建 了 近 60 a 湖泊 面积 .水 位 和 水 量变 化 时 间 序 
列 ,结合 长 时 序 气象 ,水文 观测 数据 和 土地 利用 数 
据 等 ,分 析 了 降水 .和 人 湖 径 流量 .蒸发 及 地 下 水 等 水 
量 平 衡 分 量 的 变化 特征 ,并 定量 解析 了 气候 变化 与 
人 类 活动 多 因素 对 伊 塞 克 湖 水 量变 化 的 影响 。 主 
要 研究 结论 如 下 : 

(1) 1960 年 以 来 , 伊 塞 克 湖 面积 \ 水 位 与 水 量变 
化 的 同步 性 较 强 , 呈 先 下 降 后 上 升 的 趋势 ;其 中 
1960 一 1986 年 呈 快 速 减 小 趋势 ,1987 年 起 减 小 的 趋 
势 趋 缓 , 至 1998 年 后 随 着 降水 量 增加 湖泊 水 位 和 面 
只 变化 转 为 上 升 趋 势 ;1998 年 为 伊 塞 克 湖 水 量变 化 
的 时 间 拐 点 。 

(2) 不 同时 期 水 量 平衡 及 各 要 素 的 变化 差异 显 
著 。1960 一 1986 年 伊 塞 克 湖 流域 降水 量 总 体 处 于 
较 低 水 平 ,入 湖 流 量 主 要 受 农 业 灌溉 引水 影响 持续 
减少 ;1986 年 伊 塞 克 湖 流域 降水 较 1985 年 显著 增 
加 , 且 此 后 在 波动 增加 中 保持 较 高 的 水 量 ,同时 持 
续 升 温 加 速 了 冰雪 消融 ,导致 人 湖 流量 激增 ,1986 
年 和 人 湖 流量 发 生 了 由 减少 到 增加 的 突变 。1986 年 
以 前 大 部 分 年 份 湖泊 水 量 呈 负 平 衡 ,地 下 水 持续 补 
给 湖泊 ,其 中 1976 年 地 下 水 补给 最 多 ;1986 年 起 湖 
泊 的 水 量 收 支 亏损 逐渐 减 小 ,1998 年 后 随 着 人 湖 径 
流量 和 湖区 降水 量 的 增加 , 伊 塞 克 湖 水 量 平衡 转 为 
正 增长 。 

(3) 湖区 尺度 上 ,入 湖 径 流 、 降 水 、 莹 散发 等 水 


量 平衡 分 量 的 互动 关系 决定 了 伊 塞 克 湖 水 量 的 变 
化 ;流域 尺度 上 ,气候 变化 和 人 类 活动 通过 改变 入 
湖 径流 间接 驱动 湖泊 水 量 的 变化 。 以 灌溉 引水 为 
主 的 人 类 活动 是 1960 一 1986 年 湖泊 水 量 迅 速 减少 
的 主导 因素 ,1987 年 以 来 ,气候 变化 因子 对 伊 塞 殉 
湖水 量变 化 的 影响 增强 ,降水 量 显著 增加 和 持续 升 
温 导 致 的 冰雪 消融 使 得 入 湖 流量 激增 ,是 这 阶段 湖 
泊 水 量变 化 趋势 扭转 并 持续 增加 的 主要 原因 ;2010 
年 后 湖区 蒸发 耗 散 对 湖泊 水 量 的 负 效 应 明显 超过 
入 湖 径流 和 降水 的 正 效应 ,导致 了 近期 伊 塞 克 湖 面 
吕 水 位 和 水 量 的 小 幅 波 动 减 小 。 
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Abstract: Based on the satellite data, the water level and area information of Lake Issyk-Kul was extracted, and 
the water volume was reconstruct; combined with CRU meteorological data from 1960-2020, the temperature 
and precipitation observation data from 1960-2000 at the Cholpon- Ata meteorological station and the water 
volume observation data into the lake, the lake water balance model was established to analyze the changing 
characteristics of each element of the water balance and to explore its influencing factors. The results indicated 
that: (1) Since 1960, the water volume of Issyk-Kul Lake has undergone a process of continuous decrease and 
fluctuating increase, with 1998 being the inflection point of the change; from the 1960s to the mid—1980s, the 
water volume into the lake decreased continuously mainly due to the influence of irrigation diversions, and then 
turned to an increasing trend after 1986 with the decrease of irrigation water and the increase of precipitation and 
glacial meltwater; precipitation in the lake area increased at a rate of 9.1 mm- (10a) ', and the evapotranspiration 
tends to increase significantly with increasing temperature and lake area. (2) Before the mid-1980s, Issyk-Kul had 
a negative water balance in most years, and groundwater continued to recharge the lake; since 1986, the water 
balance deficit of the lake gradually decreased, and since 1998, the positive balance has been dominated. (3) The 
interaction of water balance components such as runoff, precipitation and evaporation determines the changes in 
lake water volume, while climate change in the flow-producing areas and irrigation diversions in irrigation areas 
indirectly drive the changes in lake water volume by changing runoff; from 1960 to 1986, human activities, 
mainly irrigation diversions, were the dominant factor driving the changes in Issyk- Kul water volume, with a 
contribution of 71.6%, and since 1987 the cumulative contribution of climate change factors to changes in lake 
water volume exceeds 80%. 
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